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ABSTRACT 

We present 3D-HST, a near-infrared spectroscopic Treasury program with the Hubble Space Telescope for 
studying the physical processes that shape galaxies in the distant universe. 3D-HST provides rest- frame optical 
spectra for a sample of ~7000 galaxies at 1 < z < 3.5, the epoch when ~60% of all star formation took place, 
the number density of quasars peaked, the first galaxies stopped forming stars, and the structural regularity 
that we see in galaxies today must have emerged. 3D-HST will cover three-quarters (625 arcmin^) of the 
CANDELS Treasury survey area with two orbits of primary WEC3/G141 grism coverage and two to four 
orbits with the ACS/G800L grism in parallel. In the IR these exposure times yield a continuum signal-to-noise 
of ~5 per resolution element at Huo ^23.1 and a 5a emission line sensitivity of ^ 5 x 10~^^ erg s~^ cm~^ for 
typical objects, improving by a factor of ~ 2 for compact sources in images with low sky background levels. 
The WEC3/G141 spectra provide continuous wavelength coverage from 1.1 to 1.6 /im at a spatial resolution of 
~ Of' 13, which, combined with their depth, makes them a unique resource for studying galaxy evolution. We 
present an overview of the preliminary reduction and analysis of the grism observations, including emission 
line and redshift measurements from combined fits to the extracted grism spectra and photometry from ancillary 
multi-wavelength catalogs. The present analysis yields redshift estimates with a precision of a(z) = 0.0034(1 + 
z), or cr(v) ^ 1000 kms~^ We illustrate how the generalized nature of the survey yields near-infrared spectra 
of remarkable quality for many different types of objects, including a quasar at z = 4.7, quiescent galaxies at 
z ~ 2, and the most distant T-type brown dwarf star known. The combination of the CANDELS and 3D-HST 
surveys will provide the definitive imaging and spectroscopic dataset for studies of the 1 < z < 3.5 universe 
until the launch of the James Webb Space Telescope. 
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1. INTRODUCTION 

The investment of thousands of orbits of Hubble Space 
Telescope (HST) time has provided an incomparable imag- 
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ing legacy that has revolutionized our understanding of ob- 
servational cosmology and galaxy formation. The HST's lo- 
cation in low-Earth orbit enables high- spatial resolution free 
from the distorting effects of the atmosphere and remarkable 
sensitivity, particularly at near-infrared (IR) wavelengths, due 
to much lower background levels compared to the ground. 
As but two examples, these unique capabilities have helped 
confirm the discovery of the accelerating expansion of the 
universe (Riess et al. 2004) and helped to establish that the 
epoch at 1 < z < 3 is a critical period for the evolution of 
galaxies, when ^60% of the cosmic star formation took place 
(e.g., Hopkins & Beacom 2006; Bouwens et al. 2007) and 
the structural regularity of galaxies seen today emerged (e.g., 
Elmegreen et al. 2007; Wuyts et al. 2011) 

The interpretation of deep, high spatial resolution HST im- 
ages of galaxies at z > 1 is typically limited by the lack of 
the crucial third dimension, redshift, and other physical di- 
agnostics that can be measured from the galaxies' spectra 
such as their star formation rates (SFRs) and metallicities. 
Large spectroscopic surveys selected at optical wavelengths 
are typically limited to z ^ 1 for magnitude-limited samples 
(e.g., zCOSMOS, Lilly et al. 2007) or to color-selected, UV- 
luminous galaxies at z > 2 (e.g., Steidel et al. 1999, 2003), 
which represent only a biased subset of the M > galaxy 
population at these redshifts (van Dokkum et al. 2006). Rep- 
resentative spectroscopic samples of galaxies at z > 1 must 
be observed in the near-IR, which contains redshifted light 
emitted primarily by longer-lived stars (e.g., Maraston et al. 
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2010; Wuyts et al. 2012) as well as well-calibrated spectral 
features such as the Ha and [O 111] emission lines (e.g., Erb 
et al. 2006). However, the high sky background at these wave- 
lengths makes near-IR spectroscopic surveys expensive and 
has limited sample sizes to the order of a few dozen to ~ 100 
galaxies at z > 1.5 (e.g., Kriek et al. 2008; Forster Schreiber 
et al. 2009; Law et al. 2009; Gnerucci et al. 2011; Mancini 
et al. 2011). In order to measure the statistical properties 
of the full galaxy population at z > 1, IR- selected photomet- 
ric surveys rely on redshifts and crude spectral diagnostics 
estimated from multi-wavelength sampling of galaxy spec- 
tral energy distributions (SEDs) with broad-band or medium- 
band filters (e.g., Ilbert et al. 2010; Brammer et al. 2011). 
Current state-of-the-art photometric techniques provide a red- 
shift precision of cr/(l +z) ~ 2% (e.g., Whitaker et al. 2011) 
and stellar mass estimates with precision '^0.1 dex (Taylor 
et al. 2011), but with potentially large systematic uncertain- 
ties (e.g., Marchesini et al. 2009). These photometric mea- 
surements are insufficiently precise to study the detailed prop- 
erties of individual galaxies and their local environment: the 
photometric analyses can suffer from type-dependent system- 
atics and even the best photometric redshift measurements 
still correspond to an uncertainty of '^60 Mpc at z = 2, 
many times the correlation length of galaxies at this red- 
shift (cf. -11 Mpc; Wake et al. 2011). 

Grism spectroscopy from space provides a promising 
bridge between ground-based photometric and spectroscopic 
surveys, combining the depth and multiplexing capabilities 
of the former with the spectral diagnostics of the latter. Us- 
ing the HST-mCMOS G141 grism spanning the / and H 
bands, McCarthy et al. (1999) detected Ha emission lines at 
0.75 <z< 1.9, free from the limitations imposed by the near- 
IR atmospheric absorption and OH emission features. Yan 
et al. (1999) and Shim et al. (2009) use Ha observed with the 
NICMOS grism to estimate the SFR density of the universe 
out to z = 1.9. The line sensitivity of the early NICMOS grism 
observations is similar to that of wide-field narrow-band pho- 
tometric surveys (e.g., HIZELS z = 2.3, Geach et al. 2008), 
but the grism observations sample similar or larger cosmic 
volumes for even a limited number of pointings due to the 
fact that they cover a much broader range of redshifts. At 
optical wavelengths, deep observations with the Advanced 
Camera for Surveys (ACS) G800L grism have been used to 
confirm the redshifts of 4 < z < 7 Lyman-break and Lya- 
emitting galaxies (Pirzkal et al. 2004; Rhoads et al. 2009), 
to confirm the redshifts of passive galaxies at (z) = 1.7 from 
rest- frame UV spectral features (Daddi et al. 2005), and to 
study low-mass line emitting galaxies at < z < 2 (Xu et al. 
2007; Straughn et al. 2009) at magnitudes fainter than those 
reached by typical ground-based spectroscopic surveys. For 
large-area surveys, however, it is difficult for the ACS grism 
to compete with ground-based instruments, as much larger 
telescope apertures are available on the ground where the sky 
background at optical wavelengths is not as much of a limiting 
factor. 

The improved sensitivity and larger field size (— 2 x 2 ar- 
cmin) of the Wide-field Camera 3 (WFC3), installed on the 
HST in 2009, increases the survey speed by a factor of —20 
compared to the NICMOS grism (Atek et al. 2010) and en- 
ables a true wide-field near-IR spectroscopic survey with the 
HST that would currently be infeasible from the ground. The 
first single science pointing with the WFC3 grism, taken as 
part of the Early Release Science program (PI: O'Connell, 



GO- 11359), demonstrated the remarkable capabilities of this 
mode: van Dokkum & Brammer (2010) used the B aimer ab- 
sorption features in the spectrum of a single extremely mas- 
sive, quiescent galaxy at z = 1.9 to measure a precise redshift 
and put strong constraints on its star formation history, and 
Straughn et al. (201 1) measured redshifts for 48 emission line 
galaxies at0.2<z<2.3 down to mpoQSM = 26.9. Atek et al. 
(2010) present preliminary results of the WISP parallel sur- 
vey with the WFC3 G102 and G141 grisms, showing their 
power for detecting faint emission lines over an extended sur- 
vey area. 

In this paper we present the 3D-HST survey, an HST 
Treasury program that is providing WFC3/IR primary and 
ACSAVFC parallel imaging and grism spectroscopy over 
~625 arcmin^ of well-studied extragalactic survey fields. 3D- 
HST is designed to complement the deep, multi-epoch WFC3 
and ACS imaging of the large multi-cycle "Cosmic Assem- 
bly Near-infrared Deep Extragalactic Legacy Survey" (CAN- 
DELSy^ (Grogin et al. 2011; Koekemoer et al. 2011), by 
providing spatially-resolved rest- frame optical spectra of L > 
0.5L* galaxies over 0.5 < z < 3.5. van Dokkum et al. (2011) 
present the first science results from 3D-HST, demonstrating 
a robust and remarkable diversity within a complete sample 
of massive (> 10^^ Mq) galaxies at 1 < z < 1.5 thanks to 
the unique resolved rest- frame optical imaging and Ha spec- 
troscopy at these redshifts. In the sections below, we describe 
the observation strategy, pointing layouts, current data reduc- 
tion and analysis, and sensitivity of the survey. We present 
examples of some noteworthy spectra that illustrate the capa- 
bilities of 3D-HST and we conclude with a discussion of the 
core science goals of the survey. Magnitudes listed through- 
out are given in the AB system. 

2. OBSERVATIONS 
2.1. Filters and grisms 

The WFC3 G141 grism is the primary spectral element used 
for the 3D-HST survey. The combined transmission of the 
HST optical telescope assembly and the primary spectral or- 
der of the G141 grism ("+P^") is greater than 30% from 1.10 
to 1.65 /im, reaching a peak of nearly 50% at 1.45 jim. The 
mean dispersion of the order is 46.5 Apixel"^ {R ^ 130) 
and varies by a few percent across the field of view. The 
uncertainties of the wavelength zeropoint and dispersion of 
the G141 grism are 8A and 0.06Apixer\ respectively. For a 
full description of the calibration of the WFC3/G141 grism, 
see Kuntschner et al. (2010). Spectral features covered by 
the G141 grism include Ha at 0.7 < z < 1.5, [O I1I]A5007 at 
1.2 < z < 2.3, [0 11]A3727 at 2.0 < z < 3.4, and the Balmer / 
4000 A break at 1.8 < z < 3.1 (Figure 1). The nominal G141 
dispersion corresponds to ~1000 km s~^for Ha at z > 1; how- 
ever, in practice the resolution of the slitless grism spectra is 
determined by the physical extent of a given object (see Sec- 
tion 4.2). 

Observations with the HST grisms typically require an ac- 
companying image taken with an imaging filter to establish 
the wavelength zeropoint of the spectra (see, e.g., Kiimmel 
et al. 2009). For 3D-HST, we obtain these so-called "di- 
rect" images in the broad F140W filter that spans the gap 
between the standard / and H passbands and lies roughly in 
the center of the G141 sensitivity (Figure 1). While CAN- 
DELS will eventually provide significantly deeper imaging 

http://candels.ucolick.org/ 
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Figure 1. Throughput curves of the WFC3/G141 (primary) and ACS/G800L (parallel) grisms and the WFC3/F140W and ACS/F814W imaging filters used to 
define the wavelength reference for the grisms. The shaded bands in the top panel indicate the redshift range in which rest-frame optical spectral lines of O 11, 
H/5, O 111, and Ho; fall within the coverage of the WFC3 (red) and ACS (blue) grisms. 



of the 3D-HST fields in the F125W and F160W WFC3 fil- 
ters, the 3D-HST F140W direct images can be useful for sci- 
entific analysis in addition to calibrating the grism, as they 
reach depths competitive with even the deepest ground-based 
surveys {H ~ 26.1, 5cr) with spatial resolution ~ Of' 13 (see 
Section 3.2 and also van Dokkum et al. 2011). 

In addition to the primary WFC3 observations, 3D-HST ob- 
tains parallel ACS F814W imaging and G800L grism spec- 
troscopy. The G800L grism covers wavelengths 0.55-1.0 /im 
with a dispersion of 40 Apixel"^ and a resolution of 80 A 
for point-like sources (Kiimmel et al. 2011b). The parallel 
spectroscopy extends the Ha line sensitivity of the survey to 
z = and will provide coverage of additional lines at cer- 
tain redshift intervals where only a single line is visible in 
WFC3/G141, for example [O III] in G800L and Ha in G141 
at0.7<z< 0.8 (Figure 1). 

2.2. Survey fields 

The 3D-HST survey is allocated 248 primary+parallel or- 
bits over Cycles 18 and 19. 3D-HST will cover roughly 
75% of the area imaged by the CANDELS survey in the 
EGS/AEGIS, COSMOS, UKIDSS-UDS, and GOODS-South 
fields (Figure 2). Furthermore, G141 grism coverage of 
most of the GOODS-N field from program GO- 11 600 (PI: 
B. Weiner) is incorporated into 3D-HST as the observational 
strategy of the GOODS-N observations are nearly identical to 
that of 3D-HST. These parts of the sky are the best- studied 
extragalactic survey fields, offering a wealth of deep, multi- 
wavelength imaging and spectroscopic observations from 
large number of previous ground-based and space-based sur- 
veys. Grogin et al. (2011) provide a detailed summary of the 
ancillary datasets available in these survey fields, from X-ray 
to radio wavelengths, which are crucial for the interpretation 
of both the WFC3 imaging and grism spectroscopy datasets 
(see also Section 4.2). 

2.3. Observational strategy and mosaic layout 

The 248 3D-HST orbits are divided among 124 individual 
visits of two orbits each. The layout of the 3D-HST pointings 
is shown in Figure 2 and summarized in Table 1. In order 
to schedule 3D-HST concurrently with CANDELS observa- 
tions of the same fields over Cycles 18 and 19, initially no 
ORIENT constraints were imposed on any of the visits. Af- 
ter scheduling the ORIENTS were fixed and the positions of 
the individual pointings were optimized to provide contiguous 
mosaics and maximum overlap between the primary WFC3 



G141 and parallel ACS G800L observations, as shown in Fig- 
ure 2. Owing to this optimization, fully 90% of the G141 mo- 
saic (excluding GOODS-N) will be covered by between two 
and four orbits of the ACS grism. Within the GOODS-South 
mosaic, four additional two-orbit visits at the same orientation 
are centered on the Ultra Deep Field (UDF). The GOODS- 
South pointings outside of the area with CANDELS cover- 
age provide WFC3 grism spectroscopy of the HUDF09^^ and 
WFC3-ERS fields (see Bouwens et al. 2010). 

The first 3D-HST exposures were obtained in 2010 Octo- 
ber and the survey is nearly completed as of 2012 March: 
all pointings in the COSMOS, GOODS-S, and UDS fields 
have been observed, and two failed pointings in AEGIS will 
be re-observed by the end of 2012. The 28 grism point- 
ings in GOODS-N were completed in 2011 April.^^ Due to 
scheduling constraints, a given field will not have both com- 
plete CANDELS and 3D-HST coverage before both surveys 
are finished. For example, the two epochs of CANDELS ob- 
servations in the UDS were completed in 201 1 January, while 
the 3D-HST coverage of that field was only completed in 2012 
February. 

Each of the 3D-HST two-orbit visits with WFC3 is struc- 
tured in an identical fashion: four pairs of a short F140W 
direct image followed by a longer G141 grism exposure. The 
four pairs of direct+grism exposures are separated by small 
telescope offsets to enable the rejection of hot pixels and pix- 
els affected by cosmic-rays, as well as dithering over some 
WFC3 cosmetic defects such as the "IR-blobs" (Pirzkal et al. 
2010). The dither pattern is shown schematically in Figure 
3. The sub-pixel offsets of the dither pattern are chosen to 
improve sampling of the WFC3 PSF, which enables some re- 
covery of the image quality lost by the pixels that undersam- 
ple the PSF by a factor of 2 (see, e.g., Fruchter & Hook 2002; 
Koekemoer et al. 201 1). Due to the offset and rotation of ACS 
with respect to WFC3, the pixel subsampling of the ACS par- 
allel exposures is somewhat less than optimal (Figure 3). 

All of the WFC3/IR exposures are obtained in MULTIAC- 
CUM mode using either a SPARS50 (F140W) or SPARS 100 
(G141) read-out sequence. All of the direct exposures have 
NSAMP=5, corresponding to individual exposure times of 
203 s and a total of 812 s for each visit/pointing. Depend- 
ing on the varying usable length of an orbit, the G141 ex- 
posures have NSAMP=12-15 and the total grism exposure 

http : / /archive . st sci . edu/prepds/hudf 9/ 
The fitting analysis and specific spectra described below come from 70 
G141 pointings available as of August 1, 2011. 
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Figure 2. Layout of the 124 3D-HST pointings. Primary WFC3 F140W+G141 pointings are shown in blue with the pointing ID numbers as defined in the 
//^SrPhase II file. The locations of the parallel ACS F814W+G800L observations are shown in light green. Also indicated is the distribution of the 28 pointings 
covering the GOODS-North field from program GO- 1 1600 that are incorporated into 3D-HST. The light gray polygons indicate the footprint of the CANDELS 
WFC3 imaging, including both the "wide" and "deep" components of that survey. Note that the relative sizes of the separate fields are not shown exactly to scale. 



Table 1 

Summary of the 3D-HST survey fields 



Field 


Prog. ID 




R.A. 




Dec. 






Sky (e" s"^) 


AEGIS 


12177 


14 


: 19 : 


31 


+52 : 51 


00 


30 


0.9 


COSMOS 


12328 


10 


: 00 : 


29 


+02 : 20 


36 


28 


2.4 


GOODS-South 


12177 


03 


: 32 : 


31 


-27 : 48 


54 


32 


1.4 


HUDF09 


12177 


03 


: 32 : 


39 


-27 : 47 


01 


3b 


1.2 


UKIDSS-UDS 


12328 


02 


: 17 : 


26 


-05 : 12 


13 


28 


1.4 


GOODS-North 


11600^^ 


12 


: 36 : 


50 


+62 : 14 


07 


28 


1.0 



^ N the number of independent WFC3 pointings within the survey fields. 
^ The main HUDF09 field is covered by four visits. The flanking HUDF09-1/2 
fields have one visit each. 
Incorporated into 3D-HST; PI: B. Weiner 

times range from 4511-5111 s per visit/pointing. The ACS 
F814W direct images all have total exposure times of 480 s 
per visit, and the G800L grism images have exposure times 
ranging from 2925 to 3523 s. The ACS exposure times are 
somewhat less than those of WFC3 as a result of the larger 
overheads for reading out the ACSAVFC detectors. 

3. DATA REDUCTION 

3.1. Pre-processing 

We use as a starting point the standard calibrated data prod- 
ucts (images with the FLT extension) provided by the the 
if^Tarchive that have been processed by the calwf 3 and 
calacs reduction pipelines for the WFC3 and ACS 3D- 
HST exposures, respectively. The calibrated WFC3/IR im- 
ages have 1014x1014 pixels, with roughly 0f'128 pixer^ 
The ACSAVFC images consist of two 4096x2048 pixel ex- 
tensions with roughly Of'05 pixer^ Briefly, the calibration 
pipelines flag known bad pixels in the data quality image ex- 
tensions, subtract the bias structure determined from zeroth 
reads and overscan regions, subtract dark current, and ap- 
ply multiplicative corrections for the detector gain. Addition- 
ally, the pipelines apply a multiplicative flat-field correction to 



the WFC3/F140W and ACS/F814W direct images. The flat- 
fielding of the grism exposures is discussed below in Section 
3.2.2. Finally, bias striping (Grogin et al. 2010) and charge- 
transfer efficiency (CTE) (Anderson & Bedin 2010) correc- 
tions are applied to the ACS FLT pipeline products. A more 
comprehensive description of the WFC3 and ACS reduction 
pipelines is given by Koekemoer et al. (2011). 

3.2. Image preparation 

With the calibrated images obtained as described above, 
we perform a number of additional preparation steps in or- 
der to produce mosaics of each visit independently, which are 
composed of the four-exposure sequences described in Sec- 
tion 2.3. First we use the MultiDrizzle software (Koeke- 
moer et al. 2002) to identify any hot pixels or cosmic rays not 
flagged by the instrument calibration pipelines. We use es- 
sentially the same strategy as described by Koekemoer et al. 
(2011), though we note that we had to increase the first-pass 
detection threshold of the cosmic ray rejection to avoid re- 
jecting the central pixels of stars. The final data product after 
applying the alignment and sky-subtraction steps described 
below is an undistorted mosaic with a final scale of 0f'06 
pixel"^ (with PIXFRAC=0.8), which is justified by the opti- 
mal sampling of the WFC3 PSF by the 3D-HST dither pattern 
(Figure 3). The 5cr depth of the F140W detection images is 
if 140 ~ 26.1 for point sources within a O'.'S diameter aperture, 
varying by ±0.3 mag due to the field-dependent background 
levels (see Table 1). 

3.2.1. WCS alignment 

Small adjustments to the commanded telescope 
dither offsets are determined using the PyRAF routine, 
tweak shifts. These shifts are at most 0.1 pixels. We 
refine the WCS coordinates of the mosaic image with respect 
to a WCS reference image by matching object catalogs 
extracted from each image and fitting for shifts and rotations 
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Figure 3. Dither pattern in WFC3 pixels for each two-orbit visit, starting 
at the origin (0,0). At each position, a short F140W direct exposure is fol- 
lowed by a longer G141 exposure at the same position. The insets show the 
sub-pixel sampling of the dither pattern for the four primary WFC3 (blue) 
and four parallel ACS (green) exposures. The WFC3 pixel is fully sampled 
at half-pixel intervals, potentially allowing optimal image combination with 
interpolation. The sampling of the parallel ACS pixels (solid line) is similar, 
but not identical, to that of optimal sampling (dashed line). 

using the IRAF task, geomap. The adopted WCS reference 
images are the following //^T-based public image mosaics: 
AEGIS, ACS-/814 (Davis et al. 2007); COSMOS, ACS-Zgn 
(Koekemoer et al. 2007); GOODS-N, ACS-zsso (GiavaHsco 
et al. 2004); GOODS-South, ACS-zgso (GiavaHsco et al. 
2004); HUDF09, WFC3-ifi6o (Bouwens et al. 2010); and 
UDS, WFC3-ifi6o (Koekemoer et al. 2011). The public 
CANDELS mosaics are used when available, and future 
CANDELS data releases will be adopted for those fields 
currently using ACS data products. The derived image 
rotations are typically less than 0.1 deg, and the rms of the 
shifts matched between the catalogs is generally ~0.1 pixel 
(0f'006). Once the relative and offset shifts are determined for 
the direct images, the grism exposures are assigned the same 
shift as their preceding direct image, assuming that there was 
no shift between the two exposures (as commanded). 

3.2.2. Background subtraction and grism flat-fielding 

While the near-IR background is much lower in low-Earth 
orbit than it is from the ground, the background flux is still a 
significant component that must be subtracted from both the 
direct and grism exposures. Subtracting the background from 
the direct images is done in the following way: we subtract 
a second-order polynomial fit to each exposure after aggres- 
sively masking objects detected in the MultiDrizzle mo- 
saic and mapped back to the distorted frame using the PyRAF 
blot routine. A polynomial fit with order greater than zero 
is necessary as some exposures show structure in the back- 
ground likely caused by enhanced airglow for a particular 
pointing orientation with respect to the Earth limb. 

Subtracting the background from the grism exposures is 
more difficult than for the direct images. There is significant 
structure in the grism background that is a result of a super- 
position of multiple dispersed grism orders and the fact that 
the finite field of view of the entrance window causes differ- 
ent regions of the detector to see different combinations of the 
spectral orders. Effort has been made to produce a "master" 




Figure 4. Master background images for the G141 grism. Each master im- 
age is a masked median combination of 10 or more individual flt images 
that show similar structure in the background, such as the relative contrast of 
the dark strips on either side of the image. The F140W flat-field was divided 
out before the image combination. Black/white parts of the displayed images 
correspond to deviations of ± 5%. a) All available COSMOS pointings, b) 
GOODS-N "low", c) GOODS-N "high", and d) GOODS-N "very-high". The 
structure in the GOODS-N background images correlates roughly, but not 
exactly, with the overall background level. Using only a single master back- 
ground image, such as the one provided with the aXe calibration files, can re- 
sult in significant amount of residual structure in the background-subtracted 
images. 



grism background image created from the (masked) average 
of many grism exposures (Kiimmel et al. 2011a), which can 
then be scaled to and subtracted from a particular grism ex- 
posure. This is the approach followed by the aXe software 
package (Kiimmel et al. 2009). 

The 3D-HST fields span a wide range of celestial coordi- 
nates that result in distinct orientations of the instrument with 
respect to the zodiacal and Earth glow sources of the back- 
ground light. From the large number of 3D-HST grism obser- 
vations that have already been obtained, we find that a single 
master G141 background image is insufficient to explain the 
observed variety of structure in the background. We therefore 
correct the background in the following way. We produced 
four representative master background images that are me- 
dian (masked) combinations of subsets of 3D-HST grism ex- 
posures that have similar overall structure, determined by eye. 
These images are shown in Figure 4. The primary features in 
the grism background are the dark vertical bands spanning 
roughly 100 pixels at the left and right sides of the image. 
The relative intensity of these bands and the sharpness of their 
edge is correlated with the overall background level, but they 
can also vary at a given background level for images in differ- 
ent survey fields. There are dark horizontal bands in some of 
the master images that look like negative grism spectra. These 
features are the result of the grism "dispersing" the IR-blob 
features of decreased sensitivity, which are not physically lo- 
cated on the detector itself but rather on the Channel Select 
Mechanism in the WFC3/IR optical path (Pirzkal et al. 2010). 

The master background images shown in Figure 4 look 
very different from the master background image produced 
by Kiimmel et al. (2011a) and distributed with the G141 cal- 
ibration files. For the images in Figure 4 we have divided 
out the F140W imaging flat-field before creating the image 
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Figure 5. Flat-field correction and background subtraction for a single 3D- 
HST pointing (AEGIS- 1 1). a) A combination of the raw grism images shows 
the combined features of the background and the WFC3 flat-field features 
such as the "wagon wheel" at lower right, b) Dividing by the F140W imag- 
ing flat removes the flat-field features. The wavelength-dependent compo- 
nent of the flat-field is generally less than 2% (Appendix A), and applying a 
single flat-field greatly simplifies the data analysis, c) Subtracting by the best 
(scaled) master background image (see Figure 4) removes most of the remain- 
ing structure in the background, d) The low-level structure that remains after 
subtracting the master background is largely parallel to the image columns 
and is removed by subtracting the average along columns after all object flux 
has been aggressively masked. The top panels show the image counts aver- 
aged along columns with the full vertical plot range indicated by the numbers 
in brackets (in e~/s). For reference, a compact object with H140 = 23 has a 
peak flux of ~ 0.05 e~/s. 

combinations in order to separate the pixel-to-pixel variations 
of the flat-field from the more smoothly- varying background. 
We find that the wavelength dependence of the flat-field is less 
than ~1% across the G141 sensitivity (Appendix A). This is 
similar to the uncertainties of the overall G141 flux calibra- 
tion, and we therefore opt for a simplified approach by divid- 
ing out the single F140W flat-field from all of the G141 grism 
exposures. We then subtract the scaled master background 
image whose structure best matches that of a particular expo- 
sure, as determined from a test. Note that it is possible that 
two subsequent grism exposures within a given visit require 
different master background images depending on the varying 
background within an orbit (Appendix B). Finally we subtract 
the median level of the resulting (masked) background pixels 
averaged along image columns. While adopting the variable 
background images greatly improves the overall background 
subtraction compared to using the single image, we find that 
this last step is necessary to remove any residual structure in 
the background, which tends to be parallel to image columns. 
The complete process of first dividing by the flat field, then 
subtracting the master background image, and finally remov- 
ing the x-dependent residuals is demonstrated in Figure 5. The 
background in the final corrected grism images is typically flat 
to better than 1%. 

We have described the flat-fielding and background sub- 
traction in considerable detail as these steps are crucial for 



minimizing systematic effects in the spectra for this blank- 
field spectroscopic survey. Local background subtraction is 
not feasible for the relatively deep 3D-HST grism exposures 
because it is generally impossible to identify pixels adjacent 
to a given object that will be entirely free of flux contributed 
by nearby overlapping spectra. Small uncorrected errors in 
the background subtraction can later be manifested as con- 
tinuum break features in the extracted spectra, particularly at 
fainter magnitudes where 3D-HST can provide truly unique 
near-IR spectra currently unobtainable from the ground. 

3.3. Extraction of the grism spectra 

We use the aXe software package (Kiimmel et al. 2009) to 
extract the grism spectra in much the same way as described 
by van Dokkum et al. (2010) and Atek et al. (2010). The pri- 
mary inputs to the extraction software are a detection image 
mosaic (F140W or F814W) and the individual background- 
subtracted grism exposures generated as described in Section 
3.2. An object catalog is generated from the detection image 
with the SExtractor software (Berlin & Arnouts 1996), 
which also produces a segmentation map that indicates which 
pixels in the direct image are assigned to each object. For a 
given pixel within a particular object's segmentation map, the 
calibration of the HST grisms determines where the dispersed 
light from that pixel will fall on the grism exposures, with the 
pixel in the direct image defining the wavelength zeropoint 
for the spectrum. Thus, the grism spectrum of an object is the 
sum of all of the dispersed pixels within that object's segmen- 
tation image, or rather, the spectrum is a superposition of the 
object profile at different wavelengths offset by the grism dis- 
persion. The effective spectral resolution is therefore a com- 
bination of the grism dispersion and the object profile in the 
sense that the effective resolution decreases with increasing 
object size in the dispersion direction (see also Section 4.2 
and Figure 10). 

Because no slit mask is used and the length of the dispersed 
spectra is larger than the average separation of galaxies down 
to the detection limit of the 3D-HST survey, the spectra of 
nearby objects can overlap. This "contamination" of an ob- 
ject's spectrum by flux from its neighbors must be carefully 
accounted for in the analysis of the grism spectra. We use 
aXe to produce a full quantitative model of the grism expo- 
sures using the information in the direct image as described 
above. This is the aXe "fluxcube" contamination model, 
which makes use of the spatial information contained in the 
high-resolution HST images to model the two-dimensional 
(2D) grism spectrum. In order to generate a model spectrum 
based on the direct image, we take the observed F140W flux 
and full spatial profile within a given object's segmentation 
map and assume a constant profile and a flat spectrum in units 
of fx. As additional HST photometric bands become avail- 
able, for example the CANDELS F125W and F160W imag- 
ing covering the 3D-HST survey fields, they will be added to 
the aXe fluxcube to incorporate the wavelength dependence 
of both the flux (i.e., the color) and spatial profile into the 
model. 

The relationship between the direct image and the grism 
spectra and a demonstration of the fluxcube spectral model 
are shown in Figure 6. The images are oriented as in the in- 
dividual FLT exposures, with the grism spectra offset in the 
positive X direction with respect to the direct image and with 
wavelength increasing towards the right. The main horizon- 
tal features in Figure 6b are the "-kP^" spectral order, which 
has the greatest sensitivity (Kuntschner et al. 2010). Com- 



3D-HST grism survey 



7 



a) Direct F140W 



c) aXe Model 



b) Grism G141 



d) Model- subtracted grism 



A 

* 



A Emission 
A 0th order 



A A 



Figure 6. Model of the grism spectra based on the observed direct images and computed with aXe. Panels a and b show 50^^ x 28^' cutouts of the F140W and 
G141 observations within the GOODS-South field, with wavelength increasing towards the right on the grism panel. Panel c shows the aXe "fluxcube" model 
of the grism spectra, where the spatial profile and intensity of the spectra are determined from the direct image(s). Along with the F140W imaging, the model for 
the example shown includes color information from the F125W and F160W CANDELS imaging of this field. Panel d shows the model-subtracted grism image, 
with the image stretch increased by a factor of two compared to panels b and c. The model is quite a good representation of the data, despite the fact that no fit 
has been done; the model inputs come from the direct image(s) and the grism calibration alone. Compact features in the model are zeroth order spectra. While 
the zeroth order spectra are not perfectly modeled and subtracted (open triangles in panel d), they can generally be identified and distinguished from emission 
lines. The majority of the residuals in panel d are emission lines, indicated by filled blue triangles. Even this small cutout shows the diversity of the emission 
lines found within the 3D-HST survey, which is fully ~1600 times larger than the area shown here. 



pact, point-like features seen in both the observed and model 
images are 0* order spectra. The residuals of the model sub- 
tracted from the grism image are shown in Figure 6d. The 
model is generally a reliable quantitative representation of the 
data. It bears mention that for the present reduction there has 
been no fit to optimize the model — the relatively low level 
of the residuals demonstrates the quality and stability of the 
G141 grism calibration. Most of the compact features in the 
residuals are in fact emission lines, which are not included 
in the model. The 0*-order spectra do not always subtract 
completely, but their presence in the model allows them to be 
distinguished from emission lines. In future versions of the 
reduction, we will implement an iterative scheme to refine the 
spectral model based on the observed spectra. 

We use aXeDr i z z le to combine the four grism exposures 
(in the original distorted FLT frame) of each visit/pointing and 
extract a 2D spectrum for each object with perpendicular spa- 
tial and dispersion axes. We adopt output 0f'06 x 22 A spec- 
tral pixels, which are roughly square with respect to the driz- 
zled pixels in the direct image mosaics. These 2D spectra with 
HST spatial resolution (- 0'/l3 for WFC3/G141) are one of 
the truly unique products of the 3D-HST survey. As there 
is no slit defining a spatial axis, an "effective slit" running 
roughly parallel to the major axis of each object is not gener- 
ally parallel to the y pixel direction in the undistorted frame. 
The aXe software includes an option to account for the ori- 
entation of this effective slit to optimize the wavelength reso- 
lution of the extracted spectra (see Kiimmel et al. 2009); how- 
ever, we simply adopt the spatial axis as perpendicular to the 
spectral trace and account for the orientation and shape of the 
object profiles in post-processing analysis (see Section 4.2). 
Finally, we use aXe to extract optimally- weighted ID spec- 
tra from the drizzled 2D spectra, where the relative weights 



are determined from the object profile perpendicular to the 
dispersion axis. 

3.4. Grism sensitivity 

We have described above how the non-negligible back- 
ground emission must be subtracted from the grism expo- 
sures. In fact, this background flux from a combination of 
zodiacal light. Earth glow, and low-level thermal emission is 
the limiting factor for the sensitivity of the 3D-HST grism ex- 
posures: the read noise of the WFC3/IR detector is ~20 e~, 
while the number of background electrons per pixel in a typi- 
cal grism exposure is 1.4 e~ s~^ x 1300 s = 1820 e~. The av- 
erage background level varies for the different survey fields, 
as shown in Table 1 , depending on the field location relative 
to the ecliptic plane and the date of observation. In general, 
the observed background levels are consistent with those pre- 
dicted by the WFC3 Exposure Time Calculator (ETC). The 
angle of the bright Earth limb can vary within an orbit, and 
the two grism exposures taken within an orbit can have back- 
ground levels that differ by as much as 50% (Appendix B) and 
also different background structure (Section 3.2.2). 

We evaluate the effective continuum and emission line sen- 
sitivities of 3D-HST using a suite of simulations that is tied 
closely to the observed F140W direct and G141 grism expo- 
sures. We use a custom developed software package mod- 
eled closely after the aXeSIM package (Kiimmel et al. 2009) 
to generate a 2D model spectrum based on (1) the spatial 
distribution of flux as determined in the F140W direct im- 
age, and (2) an assumed input (ID) spectrum, normalized to 
the F140W flux. The 2D grism spectrum is then determined 
uniquely by the grism configuration files provided by STScI 
that specify how the flux in a given pixel of the direct image is 
dispersed into the spectrum of the grism image. These scripts 
will be described in more detail and released to the commu- 
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Figure 7. Simulations of G141 spectra to evaluate the 3D-HST continuum (top panels) and emission line (bottom panels) sensitivities. The input spectrum for 
each of 13 000 simulated galaxies is a flat continuum (fx) and a single (narrow) emission line at 1.3/im with equivalent width 130 A. The left panels show how the 
continuum (per 92 A resolution element) and line signal-to-noise varies with continuum magnitude and integrated line flux, respectively. The additional panels 
show how the S/N depends specifically on higher-order properties of individual galaxies: half-light radius, background level, and morphological concentration, 
defined as the ratio of the radii containing 50% and 90% of the flux in the direct image (SExtractor flux_radius; the indicated range of C corresponds roughly 
to Sersic profiles with n = 4 and n=l from left to right). In all of the panels, the S/N dependence on the properties other than the one plotted on the ordinal axis 
have been divided out, i.e., cuts in a 5-dimensional plane, S/N = f(S, R^o, Bkg., C). The shaded gray bands in the second panels indicate R50 typical of point 
sources for 0'/06 pixels. The solid red lines in the center two panels indicate a dependence of (5S/N oc 1/ v^. 



nity in a subsequent publication; for a simple spectral model 
of a flat continuum, they produce nearly identical results to 
aXeSIM and the aXe fluxcube. 

The primary advantage of the custom software is that we 
can easily modify the full input spectrum used to generate the 
model: here we assume a simple continuum, flat in units of 
fx, combined with a single emission line at 1.3/im (i.e.. Ha 
at z = 1). The emission line has a fixed equivalent width of 
(arbitrarily) 130 A, observed frame, and the overall normal- 
ization of the spectrum (and thus the integrated line flux) is 
set to the flux_AUTO flux measured by SExtractor on the 
F140W image. The result is very much like the aXe "flux- 
cube" model shown in Figure 6, however, each modeled spec- 
trum has the same line+continuum shape. We add realistic 
noise to the simulation using the WFC3/IR noise model in the 
error extension of the FLT images, which includes terms for 
the poisson error of the source counts and the read noise and 
dark current of the WFC3 IR detector. Thus, the simulations 
fully account for noise variations as a function of background 
level across all of the available pointings, and, most impor- 
tantly, for the true distribution of source morphologies as a 
function of brightness within the 3D-HST survey. There are 
approximately 13 000 objects in the simulation. 

After computing the full grism image models, we extract in- 
dividual spectra with the standard optimal extraction weight- 
ing (Home 1986) and measure the median continuum signal- 
to-noise (S/N) between 1.4-1.6/im, averaged over a typical 
G141 resolution element of 92 A. The emission line strengths 

The simulations do not, however, account for the effects of drizzling, 
which tends to smooth out some of the apparent noise because the pixel errors 
tend to be correlated (e.g., Casertano et al. 2000). 



are extracted using the technique described in Section 4.2 and 
uncertainties on the line fluxes are determined with an MCMC 
fit of the line + continuum template combination. Thus, the 
continuum and line S/N is measured in a very similar way 
as is done in the analysis of the observed spectra. Because 
of the optimal weighting in the spectral extraction, the effec- 
tive "aperture" of the extraction is likely somewhat larger (and 
variable) compared to the 1 x 3 pixel extraction window used 
by the WFC3 ETC. 

The result of these simulations is shown in Figure 7. The 
continuum and emission-line S/N depend, clearly, on the sig- 
nal itself (S), the morphological properties of the objects and 
the noise properties of the grism exposures. We explore the 
dependence of the grism S/N on the signal itself, the ob- 
ject half-light radius, R50 (SExtractor 50% FLUX_RADIUS), 
the background count rate, and a concentration index de- 
fined as C = Rso/Rgo- The dependence of S/N on all pa- 
rameters but the one shown on the ordinal axis is removed 
in each panel to isolate the contribution of each individu- 
ally. After the strength of the signal itself, the S/N de- 
pends strongly on the object size, and background level, 
with only a weak, if any, dependence on concentration. The 
brightness where the average curves shown in black in the 
right panels cross the S/N=5 threshold — mi4o = 23.1, /a = 
5.5 X 10"^'^ erg s~^ cm~^ — represents an average, not the ab- 
solute, sensitivity of the survey. The properties of the "av- 
erage" galaxy are indicated by where the curves cross S = I 
in the right three sets of panels in Figure 7: (R^q, back- 
ground, C) ~ (Of'36, 1.4 e~ s~\ 0.44). For the optimal case 
of point sources and the minimum background, one can read 
the S/N scaling directly off the figure panels: the 5a contin- 
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Figure 8. Emission line equivalent width (observed frame) as a function of 
line flux and broad-band F140W magnitude, calculated analytically assuming 
a spectrum composed of a single emission line and a continuum flat in units 
of fx . The analytic function indicated in the inset diverges as the denomina- 
tor goes to zero, that is, zero continuum and the broad-band magnitude deter- 
mined by the line flux alone. Curves are provided for the indicated line fluxes, 
which have units of erg s~^ cm~^. The high equivalent-width ([O III]A5007) 
sample of van der Wei et al. (201 1) is shown in the red points, where the mag- 
nitude was measured in the F125W filter. The equivalent width sensitivity for 
F125W (red curve) shifts down by 24% with respect to F140W, roughly the 
difference in the filter bandwidths. The equivalent width sensitivities shown 
here are uncorrected for underlying absorption, which can be of order 4 A for 
the H Balmer lines (Savaglio et al. 2005). 

uum and line limits are approximately 1.4x1.2=1.7 and 
2 X 1.2 = 2.4 times fainter, respectively, or 23.7 mag and 
2.3 X 10~^^ erg s~^ cm~^. The sensitivity limits in the high- 
background COSMOS pointings (>2 e~ s"^ Table 1 and Ap- 
pendix B) are brighter than the average limits by a factor of 
-1/0.8(0.24 mag). 

The second column of panels in Figure 7 shows that the 
continuum S/N is roughly cx l/^/R5o while the S/N of the 
emission lines goes more like oc I/R50. This effect is ex- 
pected in the sense that the continuum sensitivity only de- 
pends on the object extent in the spatial direction as any extent 
along the spectral axis is averaged out by the dispersion of the 
grism. In the case of emission lines, however, it is the pixel 
area of the emission line extent that determines their S/N as 
more extended lines are spread out over more noise from the 
background in both dimensions. Both the line and continuum 
S/N roughly follow the square root of the background counts 
(third column. Figure 7). The result that the observed trends 
are slightly flatter than 1 / ^/x, where appropriate, likely re- 
flects additional higher-order effects that affect the S/N on the 
level of —5%, whose full characterization is beyond the scope 
of the present work. 

Given the limiting line flux determined as in Figure 7, one 
can calculate analytically the limiting line equivalent width 
as a function of observed broad-band magnitude, where the 
broad-band passband contains the line and assumptions are 
made about the continuum shape. Figure 8 shows this limiting 
equivalent width as a function of magnitude in the F140W 
filter computed for different integrated line fluxes assuming a 
flat continuum (fx). We note that the input equivalent width 
assumed in the simulations described above does not figure 
into this calculation. At the faint end, muo > 24, 3D-HST is 




MAG_AUTO (F140W - H) 



Figure 9. Cumulative number counts, per WFC3 pointing. The solid black 
line shows all objects detected in the direct images in 70 analyzed point- 
ings from COSMOS, AEGIS, and GOODS -North/South fields. The blue line 
shows the 3D-HST objects matched within 1 arcsec to an object in the an- 
cillary photometric catalogs. The green line shows the objects whose flux 
at 1.4/im is contaminated at a level less than 10% by flux from their neigh- 
bors. The purple line shows objects identified as point sources based on their 
SExtractor sizes, which are clearly separated from extended sources down to 
//i40 = 23. The orange and red lines show the number counts of galaxies with 
estimated redshifts z > 1 and z> 2, respectively, from the redshift fits de- 
scribed in Section 4.2. The shaded regions show Poisson-like la confidence 
intervals computed following Gehrels (1986). 

sensitive to (observed frame) equivalent widths > 500 A. At 
the bright end, the equivalent width limit will reach a floor 
when the poisson noise of the continuum flux is similar to the 
background noise, which occurs only at muo ^19. 

3.5. Number counts 

Figure 9 shows the number of objects brighter than a given 
magnitude found in each 2x2 arcmin WFC3 pointing of the 
3D-HST survey. The number counts shown are taken from 
70 pointings taken as of 2011 August, and the number of a 
particular type of object expected in the full 149 pointings 
of the 3D-HST survey (including the GOODS-N pointings) 
is shown on the right axis. The large majority of 3D-HST 
objects are matched within If'O to objects in the ancillary 
photometric catalogs (>90% at H < 23), though the frac- 
tion of matched objects begins to decrease significantly at 
fainter magnitudes where the shallower ancillary catalogs are 
incomplete. The fraction of objects whose spectra are signif- 
icantly contaminated by flux from their neighbors increases 
with magnitude because the object surface density increases 
and faint galaxies are more likely to be close to other faint 
galaxies and also because even the extended envelopes of 
nearby brighter galaxies can have fluxes similar to faint galax- 
ies. However, the fraction of largely uncontaminated objects 
is still ~50% even at the faint limits of the survey. Objects 
with radii less than ~ Of' 18 (SExtractor 50% FLUX_RADIUS) 
are clearly separated from resolved objects down to Huq ^ 23, 
and these point sources (i.e., stars) represent the majority of 
the brightest objects. The cumulative number counts of galax- 
ies with estimated redshifts z> I and z>2 (Section 4.2) are 
also indicated in Figure 9. Based on the available pointings, 
some 7000 galaxies with z > 1 and Huq < 23.8 are expected 
within the full survey. For most magnitudes fainter than this 
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Figure 10. Demonstration of the 3D-HST redshift fitting from two objects from the COSMOS field at z > 1. The top left panels show the full SED used in 
the fit, which includes the U — 8/im photometry and the F140W grism spectra. The best-fit (unconvolved) template is shown in light blue. The top right panels 
show the redshift probability distributions, where the light gray is for the photometry alone and the darker curve is for the fit including both the spectrum and 
the photometry. The grism spectra greatly improve the redshift constraints compared to the photometry alone. The best-fit redshift of the object at right differs 
from the independent spectroscopic measurement taken from Lilly et al. (2007) (green line) by only 0.0025 x (1 +z). The bottom panels demonstrate how the 
spectral resolution of the slitless spectra is determined by the object profile. The direct images are shown in the small inset panels and the 2D grism spectra are 
adjacent to the right. The bottom panels show continuum-subtracted spectra along with the convolved and unconvolved line fits shown in red. For the object 
on the right in particular, the 2D spatial distribution of the line emission closely follows that of the direct image, which is the sum of the line -i- continuum 
across the filter bandpass. The convolved template fits are very good representations of the observed Ha-i-N II line profiles. The integrated line fluxes 2.8 and 
4.8x 10~^^ erg s~^ cm~^ for the objects at left and right, respectively, and their equivalent widths are 43 ± 4 and 91 ± 3 A. Note that while the line is clearly 
detected by the fit for the object at left, it would be difficult to pick out the line emission "by eye" from the 2D spectrum, particularly in a blind search. 



practical continuum limit, the redshift estimates will converge 
to the limit of pure photometric redshifts limited by the depth 
of the ancillary photometry. However, 3D-HST will provide 
precise redshifts and line fluxes for a smaller sample of galax- 
ies with extreme equivalent widths, such as the star-bursting 
dwarf galaxies at z > 1.5 discovered recently by van der Wei 
et al. (2011) (see also Figure 8). 

4. ANALYSIS 

The primary goal of the 3D-HST survey is to measure pre- 
cise redshifts of a relatively unbiased sample of the galaxy 
population at 1 < z < 3 to // ^ 23.8. In many cases, redshifts 
can be easily measured from emission lines, though only a 
single emission line will be observed by either the WFC3 or 
ACS grisms for most redshifts other than within a few nar- 
row redshift intervals (see Figure 1). The combination of the 
optical and IR grism spectra will provide multiple emission 
lines for many galaxies, though the sensitivity of the parallel 
ACS/G800L grism exposures is considerably less than that 
of WFC3/G141. Redshifts measured from narrow absorption 
and broad continuum features (e.g., van Dokkum & Bram- 
mer 2010) require that the systematic effects on the contin- 
uum shape be minimized, which is a particular challenge with 
the slitless grism spectra (i.e., removing the background and 
contamination of overlapping spectra). For all of the reasons 
described above, the addition of multi-wavelength photome- 
try is critical for the redshift and SED-fitting analysis of the 
grism spectra of the general galaxy population. 

4.1. Ancillary photometric catalogs 

The 3D-HST/CANDELS fields are among the best-studied 
extragalactic survey fields in the sky. All of the fields have 
extensive multiwavelength observations from X-ray to radio 
wavelengths, the details of which can be found in Grogin 



et al. (2011). For the redshift and SED fitting analysis of the 
grism spectra, we require deep optical and near-IR photomet- 
ric catalogs selected at near-IR wavelengths, corresponding 
to the rest-frame optical for the redshifts of interest. Cur- 
rently we simply match the F140W- selected objects within 
the 3D-HST fields to the deepest K-selected catalogs avail- 
able with broad multiwavelength coverage: AEGIS, COS- 
MOS, Whitaker et al. (201 1); GOODS-S, Wuyts et al. (2008); 
GOODS-N, Kajisawa et al. (2011); and UDS, Williams et al. 
(2009). These catalogs have between 8 (UDS) and 35 (COS- 
MOS) photometric bands from U through the Spitzer-IRAC 
channels at 3-8/im. The image quality of the ground-based 
catalog detection bands (Of'7-lf'O) is not well matched to that 
of the 3D-HST imaging (Of' 13). As a result, multiple 3D- 
HST objects can be matched to single objects in the photo- 
metric catalogs (cf. ~15% at H=23), which, among other 
effects, can produce artificial 3D-HST pairs in the current 
analysis. Furthermore, even the relatively short F140W expo- 
sures are somewhat deeper than the ground-based ^-band cat- 
alog detection images (e.g., ^ < 22.8 for the Whitaker et al. 
2011 COSMOS catalog, see Figure 9). For these reasons, we 
are developing PSF-matched and deblended photometric cat- 
alogs, following Labbe et al. (2006), selected from the high- 
resolution, deep 3D-HST F140W and CANDELS images that 
will provide a one-to-one correspondence between objects 
with extracted spectra and their measured ancillary photome- 
try (R. Skelton, in prep). 

4.2. Redshift and emission line fitting 

We estimate the redshifts of 3D-HST galaxies using a mod- 
ified version of the Eazy code (Brammer et al. 2008). Two 
examples of the fits are shown in Figure 10. The inputs to 
the fit are the ID grism spectra, calibrated and extracted with 
aXe, combined with broad-band and medium-band photom- 
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Figure 11. Comparison of the 3D-HST redshift measurements to spectro- 
scopic redshift measurements taken from the Hterature. The 3D-HST sample 
is selected as indicated, with H^q < 24 and also including a cut on the Eazy 
redshift quality parameter Qz <0.2, above which the quality of the template 
fits clearly decreases (among other parameters of the fit, is incorporated in 
Qz, see Brammer et al. 2008). The redshift precision is (Tnmad = 0.0034(1 -i-z) 
(see Brammer et al. 2008 for the definition of the normalized median abso- 
lute deviation, or NMAD) with 4% catastrophic outliers at redshifts, z > 0.7, 
where strong emission lines and continuum features fall within the coverage 
of the G141 grism. The sources for the spectroscopic redshifts are as fol- 
lows: AEGIS, Davis et al. (2003); Steidel et al. (2003); COSMOS, Lilly et al. 

(2007) , Brusa et al. (2010); GOODS-South, compilation from Wuyts et al. 

(2008) , Balestra et al. (2010); GOODS-North, compilation from Barger et al. 
(2008), Steidel et al. (2003), Cooper et al. (2011). 

etry from the catalogs described above when there is an object 
in the catalog matched within I'/O of a 3D-HST object. We 
subtract the quantitative contamination model directly from 
the spectra (Figure 6). Finally, we compute a normalization 
and a linear tilt to scale the G141 spectra to the / and H 
broad-band and medium-band photometry, in order to mini- 
mize residuals affecting the continuum shape that remain af- 
ter the contamination correction and background subtraction. 
At a minimum, a zeroth order normalization is necessary to 
match the flux-calibrated spectrum to the total photometric 
fluxes that have been measured within an effective infinite 
aperture (see, e.g., Whitaker et al. 2011). 

In the standard operation of the Eazy photometric redshift 
code, the model fluxes of the high-resolution spectral tem- 
plates at a given redshift are computed by convolving the tem- 
plates with the filter transmission curve of each photometric 
band. We define effective "filter curves" for each grism spec- 
tral bin (of each object individually) by convolving a boxcar 
with the width of the spectral bin with the object profile. This 
profile is taken from the direct image and averaged along the 
spectral dimension. This accounts for both the possible tilt of 
the effective slit defined roughly along the object major axis, 
which is not necessarily perpendicular to the dispersion axis 
of the slitless spectra (Section 3.3), as well as higher-order 
features in the object profile (e.g., spiral arms). In principle 
the template fitting could be done on the full 2D spectra di- 
rectly, though we opt for the relative computational simplicity 
of using the extracted ID spectra for the present analysis. An 
implicit assumption of this method is that the 2D profile is 
"gray", i.e., that the profile is the same in the line as in the 
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Figure 12. Deviations of the grism redshifts from the previously-measured 
spectroscopic value as a function of magnitude (left panel) and signal-to- 
noise of the Ha emission line (right panel). The thick black lines shows 
the average deviation in bins with equal numbers of points (N = 46) and the 
shaded gray region shows the observed NMAD scatter within these bins. For 
comparison the solid orange and red lines indicate values of constant scat- 
ter of (Tnmad = 0.0025 and 0.005, respectively, with respect to the median. 
While there is little trend in the grism redshift offset or scatter as a function 
of magnitude, the scatter decreases significantly with increasing emission- 
line strength. 



continuum across the spectral range covered by the grism. In- 
vestigating the validity of this assumption is in fact one of 
the primary goals of the 3D-HST survey (see, e.g.. Nelson 
et al. 2012), made possible by the superb spatial resolution 
provided by the HST. 

The default set of the Eazy spectral templates con- 
tain emission lines added by determining an effective SFR 
for each template and then adding lines with empirically- 
calibrated line ratios (Brammer et al. 2011, see also Ilbert 
et al. 2010). In order to fit the 3D-HST spectra + photome- 
try we remove the emission lines from the galaxy templates 
and provide separate templates for emission lines of Ha, H/3, 
[O III]A4959+5007, and [O II]A3727 individually. The non- 
negative normalizations of the galaxy templates and emission 
lines are computed simultaneously by the code, providing 
an implicit measurement of the line strengths and equivalent 
widths along with the redshift fit. Also by construction, the 
equivalent widths of the hydrogen B aimer lines are corrected 
for underlying absorption. 

For all but the most extreme velocity-broadened lines, the 
line shape is determined by the object profile. For example, 
the spectral resolution of an object with (spatial) FWHM = 
0'/5 is approximately 180A, or 1500 kms~^for Ha at z = 1.3. 
Due to the low spectral resolution. Ha and N IIA6550+6584 
are not resolved and the Ha line measurements represent the 
sum of these line species. The G141 spectral resolution tends 
to be just sufficient to produce an asymmetrical profile for the 
[O III]A4959+5007 doublet, which can help in differentiating 
it from Ha assuming that the profile of the line-emitting re- 
gion is roughly symmetric (see also Atek et al. 2010; Trump 
et al. 201 1 for examples of Ha and [O III] line profiles). 

A comparison of redshifts determined from the 
grism+photometry fits to ground-based spectroscopic 
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Figure 13. Coverage of ground-based spectroscopic redshift measurements as a function of observed (i-H) color and i magnitude. All galaxies in the available 
3D-HST pointings are shown with gray symbols, and those marked in green have ground-based spectroscopic measurements taken from the references summa- 
rized in Figure 11. The dashed lines in all panels indicate a constant i = 22.5. Analogous to the "distant red galaxy" selection of Franx et al. (2003), galaxies with 
(i-H) > 2.2 (purple lines) are nearly all at z > 1. The solid blue and red histograms indicate the fraction of galaxies with ground-based spectroscopic redshifts 
bluer and redder than this limit, respectively, with the scale shown on the axes at right. 



redshift measurements taken from the Hterature is shown in 
Figure 1 1 . The precision of the grism redshifts is excellent 
(0.0035 X (1 +z) at z > 0.7) and is an order of magnitude 
better than is typically possible with high-quality broad-band 
photometry alone (cf. 0.034 x (1 Brammer et al. 2008). 
Furthermore, the redshift precision is nearly constant over 
the full redshift range of interest 1 < z < 3, whereas the 
precision of typical photometric redshifts is generally lower 
at these redshifts as relatively fewer observed photometric 
bands sample rest-frame optical wavelengths and the B aimer 
break (Brammer et al. 2011). 

The precision of the redshift measurements depends some- 
what on the strength and availability of emission lines in the 
grism spectra. This is demonstrated in Figure 12, which 
shows the deviations of the grism redshifts from the ground- 
based spectroscopic measurements as a function of F140W 
magnitude and the S/N of the Ha emission line for objects 
at 0.8 < z < 1.4 where the line lies within the coverage of 
the G141 grism. While the redshift precision is nearly con- 
stant as a function of F140W magnitude, the redshift preci- 
sion increases from only a ~ 0.0025 for galaxies with strong 
emission lines to a ~ 0.005 for galaxies without detected 
emission lines (S/N< 5). The sample shown in Figure 12 is 
dominated by objects from the extensive catalog of Barger 
et al. (2008) in the GOODS-N field (see also Figure 13), 
which has fewer ancillary photometric bands than the other 
3D-HST fields. Therefore, it is notable that the precision of 
the continuum-only grism redshifts is still significantly better 
than has been achieved for surveys using many intermediate- 
width photometric filters (e.g.. Wolf et al. 2004; Ilbert et al. 
2009; Whitaker et al. 2011). 

Of the galaxies shown in the comparison of Figure 11, 
some 4% have redshift measurements that differ by more 
than 0.1 X (1 +z) from the spectroscopic value. The causes 
of these outliers include misidentification of single emission 
lines, poor background and/or contamination subtraction of 
the G141 spectra and apparent errors of the spectroscopic red- 
shifts themselves, where a clear emission line (or lines) is 



observed in the grism spectrum and the overall SED (spec- 
trum+photometry) is well fit by the Eazy templates. In the 
current analysis, some significant systematic redshift errors 
remain at z < 0.7, where the WFC3/G141 grism does not sam- 
ple strong spectral features. As we do not currently impose 
any constraints on the emission line ratios, these errors are 
typically caused by the fitting algorithm incorrectly adopting 
emission lines with extreme equivalent widths to fit optical 
continuum features sampled by individual, non-overlapping 
broad-band filters. These effects will be removed by future 
improvements to the fitting code and, in particular, by the in- 
corporation of the ACS/G800L optical grism spectra. 

For the z > 0.7 selection shown in Figure 11, 30% of the 
full 3D-HST sample have a measured ground-based spectro- 
scopic redshift. However, this ratio varies greatly as a function 
of object magnitude and color among the survey fields de- 
pending on the selection criteria of the spectroscopic surveys 
that cover them. Figure 13 shows the distribution of galax- 
ies with measured spectroscopic redshifts as a function of H 
(mi4o) magnitude and observed (i-H) color. Spectroscopic 
surveys selected in the optical are most complete for bright, 
blue galaxies, typically at z ^ 1 . For galaxies with H <2l and 
(i-H) < 2.2, 50% in the COSMOS and AEGIS fields have a 
measured spectroscopic redshift and this fraction reaches 80% 
and nearly 100% for the deeper surveys of the GOODS-South 
and North fields. At fainter magnitudes, optically- selected 
surveys only observe the bluest galaxies and do not sample the 
significant population of redder galaxies, typically at z > 1. 
Measuring spectroscopic redshifts for this population gener- 
ally requires near-infrared spectra, which are challenging and 
time-consuming to obtain from the ground for large samples 
(e.g., Kriek et al. 2008). 3D-HST provides near-infrared spec- 
tra of all of the galaxies in Figure 1 1 , which yield both crucial 
spectroscopic coverage of galaxies missing from typical red- 
shift surveys and rest-frame optical spectra of galaxies at z ^ 1 
that may have been observed previously only at rest- frame UV 
wavelengths (see, e.g., van Dokkum et al. 2011). 
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4.3. Example G141 grism spectra 

As a blind survey of every object that falls within a given 
WFC3 (and/or ACS) grism pointing, 3D-HST produces high 
quality spectra for a broad diversity of celestial objects. Some 
examples of the NIR grism spectra, chosen not as a represen- 
tative sample but rather to demonstrate some of this diversity 
are shown in Figure 14: 

• Figures \A a) and b) show examples of line-emitting 
galaxies. The object in panel a) has multiple com- 
ponents that are all strong [O III]A4959+5007 emitters 
with large [O III]/H/3 line ratios. The spatial extent of 
the [O III] A4959+5007 is clearly visible in the 2D grism 
spectra and closely follows the spatial profile in the di- 
rect image. A typical magnitude-limited spectroscopic 
survey would likely only be able to place a slit on one 
of the two objects and their physical association would 
be unknown. There are many such line-emitting pairs 
and groups within the 3D-HST survey that will be used 
to study the merger history of the universe out to z ~ 2. 

Panel b) shows a quasar in the COSMOS field at z = 
4.656, identified from strong emission lines of Mg II 
and C II and a strong Lyman break in the ancillary op- 
tical photometry (not shown). The template fit of the 
rest-frame UV continuum is also in remarkable agree- 
ment with the G141 spectrum, showing numerous low- 
level absorption features. The galaxy is identified as an 
X-ray source by Elvis et al. (2009), but its distant red- 
shift has not otherwise been identified; it has Zphot = 0-48 
in the catalog from Ilbert et al. (2009) and is too faint 
(/ = 22.6) to be included in the zCOSMOS-bright sam- 
ple (Lilly et al. 2007). In this sense, among other things 
3D-HST is a near-IR spectroscopic survey of X-ray (or 
radio, or 24/im, etc.) selected sources within the survey 
fields. 

• Figures 14 c) and d) show NIR spectra of extremely 
massive (M > 10^^^ M©) galaxies in the COSMOS 
and GOODS -S fields. The grism spectra covering the 
B aimer/4000 A break, even with their low resolution, 
can greatly improve the constraints on the star forma- 
tion histories of massive galaxies such as these (van 
Dokkum & Brammer 2010). The remarkable galaxy 
in panel d) shows a strong 4000 A break characteristic 
of a significantly evolved stellar population. The object 
in panel d) is taken from the CANDELS grism pointing 
within the UDF (see Trump et al. 2011). The spectrum 
is somewhat deeper (^exp = 15 ks) than typical 3D-HST 
G141 spectra, but was fully reduced with the 3D-HST 
pipeline. 

• The bottom two panels of Figure 14 show spectra of 
two more nearby objects, cool brown dwarf stars with 
spectral types T5/6 and L4. The spectral types are esti- 
mated from template fits from the SpeX spectral library 
(Burgasser et al. 2010), with the two best fits shown 
plotted along with the spectra. The T5/6 dwarf of panel 
e) is the faintest (i.e., most distant) T-dwarf by almost 
three magnitudes in H (F140W)^^. The T5/6 dwarf has 

Compared to the compilation at http://www.dwarfarchives.org. 
Note also that the T dwarf has //tot = 22.9 in the more standard H filter from 
the CFHT WIRDS photometry of the AEGIS field (Bielby et al. 201 1). 



a photometric distance of 300-400 pc using the Mh re- 
lations tabulated by Vrba et al. (2004), with the broad 
range resulting from the crude estimate of its spectral 
type. Recently, three similarly faint, distant, cool brown 
dwarfs have also been found in the WISP grism survey 
by Masters et al. (2012). 

Admittedly, the objects shown in Figure 14 are somewhat 
brighter than typical galaxies found in the full 3D-HST sur- 
vey. We provide additional example spectra in Appendix C 
down to the effective sensitivity limits of the survey. Contin- 
uum breaks can indeed be detected down to the continuum 
limit (//i4o ~ 23), though the redshift precision of such faint, 
line-free, continuum objects will likely be intermediate be- 
tween that shown in Figure 12 and that of medium-band pho- 
tometric surveys such as the NEWFIRM Medium Band Sur- 
vey (see Figure 12). Robust 5cr emission line detections are 
found at considerably fainter magnitudes at line fluxes con- 
sistent with the limits described in Section 3.4. For more 
examples of 3D-HST grism spectra, see van Dokkum et al. 
(2011) for the one-dimensional spectra of a complete, physi- 
cally defined, sample of massive galaxies at 1 <z< 1.5, and 
a sample of Ha emission line morphologies is presented by 
Nelson et al. (2012). 

5. SCIENCE OB JECTIVES OF THE SURVEY 

The 3D-HST survey described above is ideally-suited for 
studying galaxy evolution in the key epoch 1 < z < 3.5, 
providing an important complement to the large multi-cycle 
CANDELS imaging program (Grogin et al. 201 1; Koekemoer 
et al. 201 1) by measuring the crucial third dimension, redshift, 
for some 10"^ galaxies (Figure 9). Grogin et al. (2011) sum- 
marize much of the science that will be enabled by a large 
//^Tnear-IR imaging program. In a discussion that is by no 
means exhaustive, we describe below some of the science 
questions that require both high resolution imaging and the 
unique spectra that currently only 3D-HST can provide. 

5.1. What causes galaxies to stop forming stars? 

In the low redshift universe many galaxies are observed to 
be quiescent, with current SFRs only ~1% of their past aver- 
age (e.g., Pasquali et al. 2006). These quiescent galaxies tend 
to be massive early-type galaxies, forming the "red sequence" 
in the color-mass distribution of galaxies. Recent work has 
shown that at z ^ 2 many massive galaxies (M > 10^^ M©) 
exhibit spectacularly high SFRs of hundreds of solar masses 
per year, whereas others were already quiescent, particularly 
those that are extremely compact for their mass (Kriek et al. 
2006; van Dokkum et al. 2008; Brammer et al. 2009). Active 
galactic nucleus (AGN) feedback is a possible mechanism to 
suppress gas cooling and star formation (e.g., Croton et al. 
2006), but direct evidence is scarce. 

Diagnostics of quiescence can be correlated with stellar 
mass, surface density (i.e., compactness), and the environ- 
ment of galaxies on Mpc scales. These diagnostics are most 
reliably identified spectroscopically, through the strength of 
the B aimer- or 4000 A break (Z)4ooo; see Figure 14c, J) and/or 
the absence of emission lines such as Ha. For a 5cr limiting 
emission line flux of 5 x 10~^^ erg s~^ cm~^ (Section 3.4), 3D- 
HST will reach Ha and [O II]A3727 SFR limits (Kennicutt 
1998) of 2.5 and 25 Mq yr"^ at z = 1 and z = 2, respectively. If 
the simultaneous presence of quiescent and star-bursting mas- 
sive galaxies at z ~ 2 is the result of their surface density or 
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Figure 14. Some of the diversity of objects within the 3D-HST survey. The template fits from the modified Eazy fits to the spectra + photometry are shown 
in the red and orange lines in panels a-d. Panel a) shows an object in the GOODS-North field with multiple line-emitting components. Two separate spectra 
are shown extracted for the bright compact component (which itself has two close sub-components) and the fainter, more diffuse tail extending to the upper right 
of the image thumbnail. Panel b) is a quasar in the COSMOS field at z = 4.656 with strong emission lines of Mg II and C II. Panels c) and d) show extremely 
massive galaxies (10^^-^ and 10^^-^ ^o) at z ~ 2 with strong continuum breaks and no visible emission lines. The inset panels show the full 0.3-8/im SEDs 
(photometry -i- spectra) and the template fit. The bottom panels e) and f) show the spectra of T- and L-type brown dwarf stars, found in the AEGIS and GOODS-N 
fields, respectively. The two best-fitting spectral templates from Burgasser et al. (2010) are plotted on top of the spectra, with the spectral types indicated. We 
emphasize that while the selection of objects shown have spectra with particularly high signal-to-noise, none of these objects are "serendipitous" detections: that 
3D-HST provides high-quality near-IR spectra of a wide variety of classes of objects is the very essence of the survey. 
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their environment, correlations should exist between the SFR 
and these parameters. 

5.2. To what extent are galaxies shaped by their 
environment? 

The morphology-density relation (Dressier 1980) states 
that "early-type" galaxies (mostly massive quiescent galax- 
ies) are relatively abundant in dense environments such as 
groups and clusters. However, at low redshift most massive 
galaxies are quiescent regardless of environment (Kauffmann 
et al. 2003; Balogh et al. 2004), so it is difficult to deter- 
mine whether the environment provides a physical mecha- 
nism that alters the galaxies (e.g., through gas stripping), or 
whether dense environments simply are the place where qui- 
escent galaxies tend to end up. There remains some tension in 
the recent literature whether galaxy-mass-driven effects (Peng 
et al. 2010) or the properties of their host dark matter halos 
(Wake et al. 2012) are the dominant factors shaping the star 
formation histories of galaxies. 

To further disentangle the roles of mass and environment, 
epochs should be considered when "massive" did not yet di- 
rectly imply "quiescent" (van Dokkum et al. 2011). The 3D- 
HST sample will be sufficiently large to determine the rela- 
tion between SFR and environment in bins of fixed mass and 
redshift. If no environmental dependence is observed at fixed 
mass, then the relations between galaxy properties and envi- 
ronment are simply a by-product of the underlying relations 
of both quantities with mass. Even excellent photometric red- 
shifts with errors fe ^ 0.04(1 +z) have a radial error on the 
comoving distance of > 150 Mpc at z = 2, which is larger 
than the distance of the Milky Way to the Coma cluster. With 
redshift errors more than an order of magnitude smaller (Fig- 
ure 11), 3D-HST will be able to provide a sensible definition 
of the environmental galaxy density on the scale of a few co- 
moving Mpc, as well as spectroscopic diagnostics of galaxy 
"quiescence". 

5.3. How did disks and bulges grow? 

The epoch 1 < z < 3 saw the wholesale transition from 
small, star-forming clumps evident in the deepest Hubble 
images to the ordered "realm of galaxies" seen today (e.g., 
Elmegreen et al. 2007; Wuyts et al. 2012). Since z ~ 1, most 
star formation has taken place in large spiral disks, but dif- 
ferent modes may have been prevalent at earlier times, such 
as disks made of star-bursting clumps that coalesce to form 
compact spheroids directly (Dekel et al. 2009). 

If bulges formed before disks, then objects with stellar 
masses 1-5x10^^ Mq and star formation concentrated on 1 
kpc scales should be seen at z ~ 3. These should evolve into 
young compact bulges surrounded by disk-like ('^5 kpc) star 
formation at z ~ 2, and then to old bulges and regular disks 
by z ~ 1 . By contrast, if bulges formed mostly from subse- 
quent merging of disks, extended star-forming disks should 
be pervasive at z ~ 3. The relative ages of the subcompo- 
nents of galaxies and the spatial extent of line-emitting star- 
forming regions can be measured from the spatially-resolved 
HST grism spectra (van Dokkum & Brammer 2010). Ground- 
based integral-field spectroscopy has demonstrated the power 
of spatially-resolved spectroscopy for galaxies at z > 1 (e.g., 
Forster Schreiber et al. 2009; Law et al. 2009), but it has been 
typically limited to more luminous, rare objects or relatively 
small samples. 

5.4. What is the role of mergers in galaxy formation? 



Although the merger-driven growth of massive galaxies is 
a common prediction of galaxy formation models, it has been 
difficult to test at higher redshift where mergers should be 
most common (e.g., Guo & White 2008). The merger rate 
can be determined from physical pair statistics, but these have 
been difficult to measure due to insufficient spatial resolution 
and contamination by chance superpositions of unassociated 
galaxies (e.g., Williams et al. 2011; Man et al. 2012). 

With its spatial resolution of Of' 13 and spectral resolution of 
Sv ^ 1000 km s~^ (Figure 11), 3D-HST can spectroscopically 
identify true physical pairs and groups down to separations 
< 5 kpc, weeding out projected galaxy pairs (see also Figure 
14a). Within the 3D-HST sample, the pair fraction can be 
measured as a function of mass and redshift, and the fractions 
can be turned into a merger rate using models (e.g., Kitzbich- 
ler & White 2008; Lotz et al. 201 1 ; WilHams et al. 201 1). The 
mass growth of galaxies due to mergers can be compared to 
the growth due to star formation, and the sizes, densities, and 
AGN content of the merger components can be determined. 
The SFRs of the spectroscopically-identified merging pairs 
can be compared to predictions of hydrodynamical simula- 
tions, such as that mergers should play a large role in driv- 
ing star formation activity and black hole accretion (e.g.. Cox 
et al. 2006). 

6. SUMMARY 

In this paper, we present the 3D-HST survey, a 248-orbit 
Treasury program to obtain low resolution (R ~ 130) slit- 
less grism spectroscopy of ~7000 galaxies at z > 1 with 
the HST grisms. Combined with the grism coverage of the 
GOODS-North field, 3D-HST will survey 625 arcmin^ of 
well-studied extragalactic survey fields (AEGIS, COSMOS, 
GOODS-N, GOODS-S, and UKIDSS-UDS), with two or- 
bits of primary WFC3/G141 coverage across the survey and 
two to four orbits with the ACS/G800L grism in parallel. 
Short "direct" images are taken in the WFC3/F140W and 
ACS/F814W filters to provide the wavelength reference for 
the grism spectra. These images are also scientifically useful 
as they reach depths comparable to deep ground-based sur- 
veys (//i4o ^ 26.1) with ~ Of' 13 spatial resolution and sample 
wavelengths between the F125W and F160W images taken 
by the CANDELS survey. 

We have developed a custom data reduction pipeline built 
around the Multidrizzle and aXe software tools, which can 
quickly and automatically reduce any typical direct + grism 
image sequence. Of particular importance to the grism data 
reduction is the subtraction of the background, which comes 
from the combination of zodiacal light. Earth glow, and low- 
level thermal emission in the IR. We demonstrate a method 
to subtract the variable structure in the background that offers 
significant improvement over the reduction with the default 
grism calibration files. The primary products of the pipeline 
are calibrated 2D and ID spectra extracted with aXe. We 
build a quantitative contamination model based on the 3D- 
HST direct images and multi-band imaging from the CAN- 
DELS survey to account for the fact that the spectra of nearby 
objects can overlap due to the lack of slits. 

The 3D-HST survey fields provide a wealth of ancillary 
multi-wavelength observations that are crucial for interpret- 
ing the grism spectra, which frequently only contain a single 
emission line, if any. We adapted the Eazy code (Bram- 
mer et al. 2008) to measure redshifts from SEDs composed 
of both the grism spectra and matched photometry spanning 
0.3-8/im. These spectrophotometric redshifts are more than 
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an order of magnitude more precise than typical broad-band 
photometric redshift estimates, with a = 0.0034(1 The 
redshift- fitting code also provides measurements of emission 
line fluxes and equivalent widths, which van Dokkum et al. 
(201 1) use to show that there is a broader diversity of Ha line 
strengths (i.e., star formation activity) among massive galax- 
ies at 1 < z < 1.5 compared to those with similar masses lo- 
cally. 

3D-HST provides a spectrum for essentially every object in 
the field (to a magnitude limit, modulo contamination from 
overlapping spectra). We demonstrate some of the diversity 
of objects found within the survey — from high-z quasars to 
brown dwarf stars — which should not be considered "con- 
taminants" or "interlopers" but are rather natural components 
of the survey. The survey is optimally designed, however, 
for the study of galaxy formation over 1 < z < 3.5. Some 
of the science objectives that require the unique combination 
of high spatial resolution, deep near-IR (and supporting opti- 
cal) spectra include disentangling the processes that regulate 
star formation in massive galaxies, evaluating the role of en- 
vironment and mergers in shaping the galaxy population, and 
resolving the growth of disks and bulges, spatially and spec- 
trally. 

All of the raw WFC3 and ACS data for 3D-HST are im- 
mediately made public in the //^Tarchive, and we plan to re- 
lease both low- and high-level data products within 18 months 
of the completion of the observations. These data products 
will include calibrated 2D and extracted ID spectra, ancil- 
lary photometric catalogs matched to the 3D-HST sample, 
and catalogs of physical parameters derived from the spec- 
tra and photometry such as redshifts, stellar masses and SFRs. 
Much of the current reduction code is made available immedi- 
ately at http : / /code . google . com/p/threedhst /, 
to which contributions and modifications from the commu- 
nity are encouraged. News updates and data releases from the 
survey will be provided on the 3D-HST webpage at http : 
//3dhst.research.yale.edu/. 
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Figure 15. Wavelength dependence of the WFC3/IR flat-field. The left panel shows the F140W imaging flat-field, divided by the pixel area map to better show 
the pixel-to-pixel variations across the entire detector. The middle panel shows the ratio of the wavelength-dependent G141 flat-field as used by aXe, calculated 
at 1.1/im and 1.6/im. The right panel shows the ratio of the pipeline imaging flat fields for the F105W and F160W filters. The color scale shown in the center 
panel is the same for the other panels. 
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APPENDIX 
A. GRISM FLAT-FIELD CORRECTION 

Each pixel in a grism image "sees" the superposition of flux at different wavelengths coming from nearby positions on the 
sky where the grism disperses the flux to fall on top of that pixel. The standard aXe data reduction tries to account for this 
by applying a wavelength-dependent flat-field correction to pixels within the dispersed spectrum of a given object. We adopt a 
more simplified treatment of the flat-field by simply dividing by the F140W imaging flat before the background subtraction and 
spectral extraction (Section 3.2.2). In doing so, we also separate multiplicative flat-field and additive background terms from the 
background subtraction. However, this technique ignores any wavelength dependence of the flat-field, which is shown in Figure 
15. 

The left panel of Figure 15 shows the F140W flat itself, with the large-scale variation resulting from the variable pixel areas^^ 
removed for display purposes. The center panel shows the estimated ratio of the flat-field at the blue to red edges of the G141 
sensitivity, computed from the wavelength-dependent flat-field images used by aXe. The right panel of Figure 15 shows the 
ratio of the pipeline imaging flats for the F105W and F160W filters. For both estimates, the wavelength dependence of the 
flat-field across the G141 sensitivity is generally less than ±1% outside of the "wagon- wheel" feature in the lower right comer 
of the detector. The color dependence of the aXe flat-field includes some high-frequency structure not seen in the ratio of the 
imaging flats, and also appears to underestimate the decreased sensitivity in the wagon wheel at blue wavelengths. Our simplified 
treatment of dividing by the single F140W flat-field can therefore result in flat-field errors of ~5% in this part of the detector, 
however, dividing by this "average" flat is generally sufficient to flatten the background and reduce systematic effects caused by 
the background subtraction (see Figure 5). 



http : //www. stsci . edu/hst /wf c3 /pam/pixel_area_maps 
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Figure 16. Background levels for all of the individual exposures in the four primary 3D-HST survey fields. The four exposures within a given visit are shown 
with circles connected by lines, showing that the background can vary significantly within a single visit (e.g., in GOODS-S). Histograms indicate the overall 
distribution of background levels in each field. 



The background level in the 3D-HST grism exposures varies significantly from 0.5 electrons/s to 3 electrons/s depending on 
the field and the relative orientation of the instrument with respect to the zodial light and Earth glow sources of background 
emission. The orientation of the bright limb can also vary within a single orbit, resulting in different overall background levels 
and also, occasionally, different 2D structure within the background (Section 3.2.2). Summarized in Table 1, we present here the 
full distribution of background levels for all of exposures in the four primary 3D-HST survey fields in Figure 16. The COSMOS 
field shows the largest variation, with a majority of the pointings having background levels in excess of 2 electrons / s. There are 
a number of GOODS-S pointings in which the background levels vary by roughly a factor of two (0.8-1.5 electrons / s) within a 
single visit. 



The objects with spectra shown in Figure 10 are shown to demonstrate the remarkable quality of spectra that can be achieved 
with the WFC3/G141 grism for moderately bright targets (though even those magnitudes are currently difficult to reach from the 
ground). However, the objects in Figure 10 are brighter than the typical object that will be found in the full 3D-HST survey (e.g.. 
Figure 9). Here we provide a gallery of spectra at continuum (Figure 17) and emission-line (Figure 18) brightnesses extending to 
the faint limits probed by the survey (Section 3.4). It bears noting that there appear to be significant narrow absorption features in 
some of the continuum spectra shown in Figure 17. While the grism is capable of resolving absorption lines (e.g., van Dokkum 
& Brammer 2010), the narrow features in the faint spectra shown in Figure 17 are more likely the result of correlated pixel noise 
that arise from the drizzling process. It is also possible that such correlated noise could cause spurious emission line detections, 
but the fact that the line shape predicted by the object morphology is in good agreement with the observed line shape in Figure 
18 demonstrates that, at least for the objects shown, even faint lines can be robustly identified. 



B. BACKGROUND VARIATIONS 



C. EXAMPLE SPECTRA OF TYPICAL OBJECTS IN 3D-HST 



3D-HST grism survey 



19 




10^ 5 X 10^ 



10000 11000 12000 13000 14000 15000 16000 17000 

A [A] 





COSMOS- 17-G141 _0045 1 
z= 1.739, mi4o = 22.1 




11000 12000 13000 14000 15000 16000 



A [A] 




10000 11000 12000 13000 14000 15000 16000 17000 10000 11000 12000 13000 14000 15000 16000 17000 



A [A] A [A] 

Figure 17. Example spectra of more typical (i.e., fainter) galaxies at z > 1.7 than those shown in Figure 10. The components of each individual object plot are 
as described in Figure 10. The galaxies shown here have well-detected continuum (B aimer/4000 A) breaks and weak, if any, emission lines. The pointing where 
each object is found is indicated in the inset labels ("MARSHALL" is a CANDELS supernova pointing in the UDS). Where emission lines of O IIIA4959+5007 
are clearly detected, the combined equivalent width and integrated flux of the doublet is indicated, in units of A and 10~^^ erg s~^ cm~^, respectively. 
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Figure 18. Example spectra of line-emitting galaxies where the line signal to noise, 3 < S/N < 8. The components of each individual object plot are as described 
in Figure 10. The equivalent widths and integrated line fluxes of Ho; and O IIIA4959-I-5007 emission lines are indicated, in units of A and 10~^^ erg s~^ cm~^, 
respectively. 



